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The growth dynamics of rigid biopolymers, consisting of N parallel protofilaments, 
is investigated theoretically using simple approximate models. In our approach, the 
structure of a polymer's growing end and lateral interactions between protofilaments 
are explicitly taken into account, and it is argued that only few conformations are 
important for biopolymer's growth. As a result, exact analytic expressions for growth 
velocity and dispersion are obtained for any number of protofilaments and arbitrary 
geometry of the growing end of the biopolymer. Our theoretical predictions are 
compared with a full description of biopolymer growth dynamics for the simplest N = 
2 model. It is found that the results from the approximate theory are approaching 
the exact ones for large lateral interactions between the protofilaments. Our theory 
is also applied to analyze the experimental data on the growth of microtubules. 
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I. INTRODUCTION 

Rigid biopolymers such as microtubules, actin filaments and intermediate filaments are 
major components of cytoskeleton and cellular environment. They play a fundamental role 
in biological systems by supporting cellular transport, cell motility and reproduction . ^i^i^i^ 
Many cellular processes essential for life are driven by polymerization/depolymerization 
dynamics of these biopolymers. Therefore, a full theoretical description of growth processes 
is clearly needed in order to understand mechanisms and principles of cell functioning. 

Microtubules are rigid, hollow tubular biopolymers made of parallel protofilaments ar- 
ranged in circular array.^^ Each protofilament is a linear polymer chain consisting of al- 
ternating a- and /3-tubulin subunits. The number of protofilaments varies between 10 and 
15 for microtubules from different species, but typically most of them have = 13 protofil- 
aments. Lattice structure of microtubules can be viewed as 3 parallel helices, the so-called 
3-start helix.i^ It also shows a discontinuity or seam, and a functional role of this lat- 
tice feature is unknown. The dynamics of microtubules features an unusual phenomenon 
of alternating between the growing and shrinking phases which is termed dynamic instabil- 
ity. Actin filaments are another example of rigid biopolymers. They can be described as 
two-stranded helices in which each actin monomer contacts four other monomers, with the 
strongest interaction along the strands.—. 

Recent experiments^"^ have provided extensive measurements of growth dynamics of 
actin filaments and microtubules under the effect of external forces. A number of the- 
oretical models aimed to describe the dynamics of growing rigid biopolymers have been 
proposed . Several studies utilized polymerization ratchet models, which as- 
sume that thermal fluctuations at the tip of growing biopolymers control the growth 
dynamics.'S*S*i2*ii In these models a rigid biopolymer is viewed as consisting of A^ inde- 
pendent and no^-interacting parallel protofilaments. A different approach is to describe the 
biopolymer's growth dynamics phenomenologically, by considering the overall polymeriza- 
tion and depolymerization processes and neglecting the microscopic details.iSiii Although 
current theoretical models provide a reasonable description of many aspects of microtubules 
and actin filaments growth, there are still many open questions.— The general deficiency 
of current theoretical approaches is the fact that they mainly ignore the microscopic struc- 
ture and geometrical properties of a biopolymer's lattice, and they also neglect the lateral 
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interactions between protofilaments. 

The purpose of this paper is to investigate the growth dynamics of rigid multi-filament 
biopolymer by taking into account a complex structure of polymer's growing end, its geomet- 
rical properties and interactions between parallel protofilaments. We show that at realistic 
conditions only few polymer configurations contribute into the overall dynamics. It allows us 
to develop a simple approximate one-layer model for rigid biopolymers growth, for which we 
obtain explicit expressions for velocity and dispersion for any number of protofilaments and 
for arbitrary geometry of the polymer end. The essential issue of this approximate treatment 
proves out to be an evaluation of analytical expressions for the asymptotic (long-time) mean 
growth velocity, V, 

^ Vt, (1) 
and for dispersion (or effective diffusion constant), D, 

[{^\t)) - mf] ■ (2) 

Here, x(t) stands for a coordinate of the biopolymer 's tip at time t that grows linearly with 
time at stationary-state limit. Dispersion D is a natural measure of fluctuations of growth 
dynamics. Note, that although in this paper we aim to describe the microtubule dynamics, it 
could also be used to analyze the growth of other rigid biopolymers, such as actin filaments. 

The paper is organized as follows: the approximate model of the growth of biopolymer 
consisting of protofilaments is presented in Sec. II, while in Sec. Ill the full description 
of growth model with N = 2 protofilaments is given and compared with the approximate 
approach. Sec. IV compares the predictions of the approximate model with the exact growth 
model for N = 2 case, and in Sec. V our theory is used to describe the real experimental 
data on the growth of microtubules. Sec. VI summarizes all results and concludes our paper. 

II. ONE-LAYER MODEL 

Consider a growing rigid polymer consisting of N protofilaments as shown in Fig. 1. The 
building block of such polymers is a monomer subunit of length d. For microtubules this 
is an a-/5-tubulin heterodimer, which has the length of 8.2 nm.i^ The protofilaments are 
parallel to each other but shifted by arbitrary distances. There are chemical interactions 
between protofilaments that bind them together in a polymer lattice. Protofilaments are 
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labeled in a such way that the seam in the polymer lattice is always between protofilaments 
1 and A^. There is an infinite number of possible polymer configurations, which can be 
described by a set of numbers, {ai,a2, . . . ,aAr}, where aj is a coordinate of the tip of 
the protofilament j. For labeling different configurations, let us choose a moving coordinate 
system where the origin is always at the tip of the leading protofilament, i.e., aj = if the 
protofilament j is the leading one: see Fig. 1. 

When a free monomer approaches the growing biopolymer molecule, it can attach to any 
of protofilaments. Define Uj (wj) as a rate of attachment (detachment) to the protofila- 
ment j. These rates are related by the following thermodynamic expression, 

— = exp[-{g^ + gh + gim)/kBT], (3) 

Wj 

where is a bond energy due to head-to-tail binding (vertical), Qh is the energy due to 
lateral (horizontal) interaction, and Qim is the free energy of immobilizing the free monomer 
into the rigid lattice:>i^ see Fig. 1. For binding to any protofilament the value of longitu- 
dinal energy is the same, however, contributions to the horizontal energy gh might be 
different because the local environment of each protofilament tip is different. For polymer- 
ization/ depolymerizat ion rates we can easily write 

Uj oc exp[-e{g^ + gh + gim)/kBT], Wj oc exp[-(6' - l){g^ + gh + gim)/kBT]; (4) 

where the coefficient 6, < 6 < 1, refiects the value of the activation barrier for the process 
of monomer binding. The exact values of the lateral interaction energy gh for microtubules 
are unknown, but there is an estimate based on computer simulations of a stochastic model 
of microtubule assembly dynamics, which gives gh approximately between -3 and -6 ksTi^ 
It is obvious that the monomer could attach much faster to some protofilaments when more 
lateral bonds are made. That leads to observation that the growing biopolymer can be 
found preferentially in some particular configurations. This is the basis for our approximate 
model. 

In our model we assume that the growing biopolymer can only be found in configurations 
with distances from all protofilament tips less than d, i.e., aj < d for all j, as shown in 
Fig. 2. It means that all protofilament ends are within one monomer distance from the 
leading protofilament. There are N such configurations, because each protofilament can be 
the leading one only once. This is a one-layer model of rigid biopolymer growth. 
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Consider the dynamics of biopolymer growth in the one-layer model. For the biopolymer 
in the configuration j, i.e. when the tip of the protofilament j is at the origin, the incoming 
monomer can only attach to one protofilament, which is the furthest away, and the polymer 
configuration transformed into another one. The whole new layer of the monomers of length 
d is added to the biopolymer only when the system goes sequentially through all allowed 
N configurations, and we return to the same configuration j. Thus, the growing biopoly- 
mer advances from the given configuration to the same one, but only shifted by distance 
d, through the sequence of states. Then the process repeats again and again. From 
mathematical point of view, this dynamic description can be mapped into the motion of a 
single particle on periodic one-dimensional lattices.— Recently, similar approach has been 
used successfully to describe the dynamics of motor proteins . 

This mapping allows us to obtain exact and explicit expressions for the stationary-state 
growth velocity and the dispersion, or effective diffusion constant, for the rigid biopolymer 
with any number of protofilaments and for any arbitrary set of shifts {ai,a2,--- ,(1^} in 
terms of rate constants {uj,Wj}. The equation for the steady mean velocity of growth is 
given byiSiiliS 

where d is the size of the monomer subunit, while 

N / N-1 k \ 

j=i ^ \ k=i i=i / 

Here we also used the periodicity of the one- layer model, i.e., Uj±N = Uj and Wj±iy = Wj. 
The expression for the dispersion is similar but more complex: ^''i^^i^^ 

d {VSN + dUN 1 

I. AT 

where the auxiliary functions are given by 

N N N 

j=l 1=1 j=l 

with supplementary coefficients 



Sj ^ 

V fc=i i=i ^j-* 
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In the simplest case of the growth of the biopolymer with N = 2 protofilaments these 
expressions are reduced to 

V = d "^"^-^^^^ , (10) 

Ui + U2 + Wi + W2 

for the mean growth velocity, while for the dispersion it yields 

j^_d'^ uiU2 + wiW2 - 2{V/df 

2 Ui+ U2+ Wi+ W2 



A. Effect of External Loads 



Forces produced by growing microtubules are crucial for understanding mechanisms of 
cellular motility and cellular transport .-^^^^ To investigate force production in experiments 
the growing microtubules are put under external loads which slow down the growth.^^ The 
experimentally used external loads are hard walls and/or optical trap systems. These studies 
provide a valuable information on biopolymer growth mechanisms and cellular motility, and 
any theoretical description should account for the effect of external forces. In the one-layer 
model this can be easily done. 

Consider a microtubule which is growing against an external force F. This force acts 
locally only on the leading protofilament. An example will be a hard wall positioned at the 
tip of the leading protofilament. When a monomer is attaching to the protofilament j, the 
microtubule produces the work and it is equal to F{d — aj). Then the rates of polymerization 
and depolymerization should be modified as follows, 

Uj{F) = Uj{0)exp[-e+F{d - aj)/kBT], 

Wj{F) = Wj{0)exp[+9jF{d - aj)/kBT] (12) 

where d — aj is the microtubule length change for monomer binding to the protofilament 
j, and 6j' and 6~ are the load- distribution factors. These factors reflect how the external 
force affects the activation energy for attachment and detachment processes of the monomer 
subunit.— Also, the load-distribution factors may have different signs, but we certainly 
expect that the overall factor 

N 

= Y^^el + 9T), (13) 

to be positive, implying an opposition to growth of the biopolymer molecule. The force- 
dependent rate constants can be substituted then in Eqs. ^ and ((Zj) to obtain the load- 
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dependent explicit expressions for the velocity and dispersion, and thus providing a full 
description of biopolymer growth under external forces (see also Fig. 6). 

If the external force opposes to polymerization, then there is a special value of force, 
termed stalling force, Fs, at which the force-dependent mean growth velocity is approaching 
zero. This is an important characteristic of real biopolymers.— In our approximate model we 
can easily calculate the stalling force, 

Fs^tf,.§l^. (14) 

Similar expressions have been used successfully for the description of motor proteins 
dynamics.—*^ 



B. Comparison with Phenomenological Description 

It is very interesting to compare our approximate approach with phenomenological mod- 
els which dominate in the field of microtubule dynamics. According to a phenomenological 
description, the growth rate is determined by the balance between polymerization and de- 
polymerization processes, 

V = -^{KnC- Kff), (15) 

where c is the concentration of free tubulin-GTP subunits, and kon and koff are average 
rate constants for polymerization and depolymerization, respectively. This simple picture 
suggests that there is a linear dependence of microtubule growth on tubulin concentration, 
at least, in the regime where the dynamic instability can be neglected. 

In the one-layer model, binding rate constants are also proportional to the concentration 
of free monomers, Uj = kjC. However, the linear dependence of the mean growth velocity 
is only valid at large concentrations of monomers, while at low concentrations the behavior 
is different. The mean growth velocity has a power law dependence on the concentration 
at this regime. It can be easily seen by analyzing Eq. ()10|) for the case of the biopolymer 
growth with N = 2 biopolymers. At large c, one can obtain V oc -^^^c^ while at low c 
the mean growth velocity has a quadratic dependence. Similar deviations are also found for 
dispersion as a function of the free monomers concentration. We expect that in general 
case the behavior is analogous. 
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The observation that the biopolymer growth velocity may deviate from the simple linear 
dependence given in Eq. (fT3j) has serious consequences since all experimental measurements 
have been based on this phenomenological picture. It means that the values of measured 
rate constants kon and koff may be different for different tubulin concentrations. This could 
explain a large concentration variability in growth rates of microtubules. ^'^^ This question 
requires a careful experimental investigation. 

III. GROWTH DYNAMICS FOR BIOPOLYMERS WITH iV = 2 

PROTOFILAMENTS 

How well the one-layer model approximates the full growth dynamics of rigid biopolymers 
is an open question. However, we can answer this question in the simplest non-trivial case of 
the biopolymer with N = 2 protofilaments, where an exact solution can be found for the full 
dynamical description. A comparison between the approximate and the complete dynami- 
cal models provides an important physical insight into the mechanism of rigid biopolymer 
growth, and allows us to understand the applicability of the one-layer model. 

When the monomer is binding or dissociating from the growing biopolymer with N = 2 
protofilaments, there are four different situations: see Fig. 3. The monomer can attach with 
the rate Ui to the protofilament 1 of the configuration {a, 0}, where a, < a < d, is the shift 
between the protofilaments, or it can dissociate with the rate Wi from the protofilament 1 of 
the configuration {0, d — a}, as shown in Fig. 3a. These rate constants are related through 
the thermodynamic expression 



where Qh is the the energy of lateral (horizontal) interaction per one subunit, is the 
longitudinal (vertical) bond energy, and Qim is the entropic term which describes the energy 
of immobilization. Then the factor 2a /d gives the fraction of the lateral bond created after 
the monomer binded to the biopolymer. 

Similarly, the monomer can attach with the rate U2 to the protofilament 2 of the con- 
figuration {0,d — a}, or it can dissociate with the rate W2 from the protofilament 2 of the 
configuration {a, 0}, as shown in Fig. 3b. The relation for these rates is given by 




(16) 




(17) 
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Rate constants U3 and W3 describe the polymerization/depolymerization processes, which 

do not change the length of the polymer, as illustrated in Fig 3c. Meanwhile, the attachment 

and detachment rates from the leading protofilament for all possible configurations are given 

by u and w, respectively, as shown in Fig. 3d. These rate constants can be described by the 

following expressions, 

U3 [ 2gh + gv + gim\ u [ gv + gim\ 

— = exp -— , and - = exp — — . (18) 

W3 \ ksT J w \ ksT J 

For convenience, let us define 7 as a new energy parameter that specifies lateral interactions, 

7 = exp {-gh/ ksT) . (19) 

Then, using Eq. (|T^ . we obtain a set of simpler relations between the rate constants, 
namely, 

ui/wi = ulw^i^"!^ 

U3/W3 = u/wj^ (20) 

Let us define P{a + id, 0) as a probability to find the system in a configuration, where 
the protofilament 2 is the leading one and the tip of the first protofilament is at distance 
a+id, {i = 0, 1, ...). Similarly, we define P(0, d—a+id) as a probability to find a configuration 
with the tip of the second protofilament at distance d — a + id, {i = 0, 1, ...) from the tip 
of the first protofilament, which is now the leading one. The overall kinetic scheme of the 
system is shown in Fig. 4. Because of the symmetry, at stationary state the overall fiux 
through the system is equal to zero. Then the following relations for probabilities are valid, 

P{a + td,0)=(^^^\ P{a,Oy, (21) 

\U3 + Wj 

P(0, d-a + id)= (^]^^^^ ^(0, d-a); (22) 
P(0,d-a) = ^^^^^P(a,Oy, (23) 

U2 + Wi 

with a normalization condition 

00 

'^[P{a + id,0) + P{0,d~a + id)] = l (24) 

2=0 
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It can be argued that the parameter = (3 is less than 1. This follows from the 

observation that u^> u and < w, i.e., the monomer binds faster and stronger to interior 
of the polymer than to the tip of the leading protofilament. After summation of geometrical 
series in Eqs. (j7T|) and (|^. we obtain 

P(a,0) = ^^ and P{0, d - a) = (25) 

where we also defined 

U1 + W2 , . 

a = . (26) 

U2 + Wi 

It is interesting to note that 1 —P{a, 0) — P(0, d — a) gives the fraction of configurations that 
are neglected in our approximate one-layer model. Using Eqs. , the simple calculation 
for this quantity yields /?. 

The mean growth velocity can be written as a sum of several terms, namely, 

V = du + {d — a)uiP{a, 0) + au2P{0, d — a) 

-dwil - P{a, 0) - P(0, d-a))-id- a)wiP(0, d-a)- aw2P{a, 0). (27) 

The first positive term corresponds to adding a monomer to the leading protofilament and 
increasing the length of the polymer by d. This may take place at any polymer configuration. 
The second and third positive terms reflect the addition to "one-layer" configurations, i.e., 
where tips of protofilaments are at the distances less than d from each other. Similarly, 
the first negative term gives the contribution from dissociation of the monomer from the 
leading protofilament. It shortens the polymer by distance d for all configurations except 
"one-layer" configurations. The last two terms represent negative contributions to the mean 
growth velocity from the "one-layer" configurations. 

Substituting expressions (^3)) into the equation ^T7\i yields the final expression for the 
mean growth velocity of the rigid biopolymer with N = 2 protofilaments, 

v^dL-wPHi-^) |. (28) 

1^ Ui + U2 + Wi + W2 ) 

Similar calculations can be performed for the dispersion, or effective diffusion constant. 

In the "one-layer" model for biopolymer growth with N = 2 protofilaments, there are 
only two configurations, (a, 0) and (0, d — a). The system can transfer from the configuration 
(a, 0) to the configuration (0, d — a) with the rate ui + W2, or it can go backward with the 
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rate U2 + wi. At t —>■ 00, the probabilities to find the system in different configurations can 
be easily calculated to give 

1 a 
P(a,0) = ^^, and P(0,d-a) = ^^. (29) 

These equations can also be obtained from the general expressions for the full dynamic 
case when (3 = 0, i.e. when we neglect all configurations except "one-layer" configurations. 
Then the mean growth velocity can be easily calculated as follows, 

V = d "^"^-^^^^ (30) 

Ui + U2 + Wi + W2 

Now, we can compare the predictions of the full dynamics and "one-layer" = 2-model 
for different values of the parameter 7, which reflects the energy of lateral interactions 
between the protofilaments. We can rewrite the thermodynamics expressions (|2m) in the 
following form, 

= U77^i-"/'^, W2 = M;7-^2-(l-a/rf)^ = W-f^'-\ (31) 

Coefficients /i, /2 and reflect the different values of activation energies for specific poly- 
merization and depolymerization processes. These coefficients may be realistically estimated 

as 

-a/d< fi<a/d, -(1 - a/rf) < /a < (1 - a/rf), -1 < /s < 1 (32) 

It means that the monomer attaches faster to the place where the stronger lateral bond is 
created. Similarly, the dissociation is slower if larger lateral bond should be broken. 

For illustration purposes only, we consider the case when /i = /2 = /a = 0. The 
calculations for other possible values of parameters /i, /2 and /a produce qualitatively 
similar results. For the model with full dynamics we obtain the following expression for the 
mean growth velocity, 

In the "one-layer" model the velocity is given by 
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To compare theoretical predictions it is convenient to analyze the ratio of velocities, which 
gives us a measure of deviations between approximate and exact approaches. 
For the ratio of growth rates in two models we obtain 



^one— layer 1 




(35) 



The ratio of velocities for two models is also plotted in Fig. 5 for different values oia/d. The 
simple analysis indicates that for any values of shift between the protofilaments the ratio of 
velocities is approaching 1 for large 7. The convergence is better with increasing values of 
a/d and reaches the maximum for a/d = 0.5, i.e., when protofilaments are shifted by a half 
of the subunit length. But even for relatively small shifts, a/d = 0.1, and realistic values 
of 7 ~ (energy of lateral interaction between protofilaments is of order 10 /c^T), the 
deviation of the approximate "one-layer" model from full dynamics description is less than 
10%. Thus the "one-layer" model provides a very good approximation for the full dynamics 
description of the growth of rigid biopolymer with N = 2 protofilaments. 



IV. APPLICATION OF THE ONE-LAYER MODEL FOR THE DESCRIPTION 
OF EXPERIMENTS ON MICROTUBULE GROWTH 

The "one-layer" model can be easily used to describe the growth of real microtubules. 
To illustrate that we apply the approximate theory to describe the experimental data of 
Dogterom and Yurke.^ In these experiments, the growing microtubules encountered a rigid 
microfabricated barrier, and the external forces have been calculated from buckling shapes 
of biopolymers for different growth velocities. 

To describe this experimental force-velocity relation we used the following parameters: 
the rate constant u for the process of association of a tubulin subunit to the tip of the leading 
protofilament with a creation of the longitudinal (vertical) bond only, the constant w for the 
dissociation process from the leading protofilament which only breaks the longitudinal bond 
(see Fig. 3d for the biopolymer with N = 2 protofilaments); the set of load distribution 
factors 9j^l < j < 13; and the free energy of creation of lateral (horizontal) bond, qh- All rate 
constants Uj and Wj can be expressed in terms of m, w, and the parameter 7 = exp(— gf/j/ ksT). 
To simplify the analysis, we also assumed that all load distribution factors are equal except 
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9i and 6'i3. It reflects the fact tlie external force may affect differently the rates of removing 
or adding a sub unit at the two adjacent protofilaments located at the seam. 

The growth velocity has been calculated using Eqs. (0) and The effect of external 
forces has been taken into account by using expressions ()12|). The rates Uj and wj have been 
expressed in terms of rates u, w, and the parameter 7 in a way similar to = 2 case. The 
resulting force-velocity curve is shown in Fig. 6. We find that the parameters 

u = 8.3s~\ w; = 355 s"\ = 1 for all j, 7 = 100, (36) 

provide a very satisfactory optimal fit to Dogterom and Yurke experimental data. Further- 
more, this fit via the relation ()14|) allows us to predict the stall force Fs — 5.5 pN. 

As was explained above, our theoretical method also provides a connection with a phe- 
nomenological description. From the experimental fit we calculated the rates kon and koff, 
which are observable rates of polymerization and depolymerization. Because the individual 
rates of association depend on concentration of free tubulin-GTP molecules, at 
the limit of very large concentrations the mean growth velocity in the one-layer model is a 
linear function. This allows us to estimate kon and koff as follows, 

N 

kon = ^/ and koff = konCe, (37) 

where Ce is the critical concentration of tubulin monomers below which the biopolymer 
growth cannot happen. It can be easily calculated from Eq. (jSJ, Cg = Y[f=i'^j/^j- For 
experimental conditions of Dogterom and Yurke^ it is given by Ce — 10.6 fiM, which is 
only slightly larger than 2-5 /iM, obtained in other experiments on microtubule assembly.— 
Calculations from the experimental fit yields the following values for phenomenological rates, 

kon = 3.0 s^VM"^ and koff = 31.8 s'K (38) 

These estimates are in a good agreement with independent experimental observations, which 
put kon in the range of 2-10 s^^/iM^\ while the experimental spread for koff is much larger, 
between 0.1 and 45 s~^. 

The "one-layer" model incorporates many microscopic properties of microtubules and it 
allows us to estimate some thermodynamic and structural properties of these biopolymers. 
The parameter 7 is associated with the lateral energy of interaction between the protofil- 
aments. From the experimental fit we conclude that this energy is given by gh — —^ksT, 
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which is in excellent agreement with the only available estimates of -(3.2-5.7) ksT, obtained 
from computer simulations of the stochastic model of microtubule assembly dynamics.—. 
Note, however, that reasonable fits could also be produced with values of 7 ranging from 10 
to 10^ 

Meanwhile the ratio of rates u/w provides information on the energy of longitudinal bond 
in microtubules. At standard conditions (1 M solution of tubulin-GTP) this energy can be 
calculated from 

gv = -keT In — - Qi^, (39) 
w 

where c = 25 fiM is the concentration of free tubulins in experiments of Dogterom and 
Yurke.=Si The parameter is a standard free energy of immobilization of a tubulin subunit 
in the polymer lattice, which was estimated to be in the range of 12-18 kBTi^. Then, using 
the fitted values of rates u and w we obtain that g^j ~ —(19 — 25) ksT. It is interesting 
to note that longitudinal interactions in microtubules are much stronger than the lateral 
interactions between the protofilaments. These estimates are in excellent agreement with 
results of computer simulations of stochastic models of microtubule growth.— 



V. SUMMARY AND CONCLUSIONS 



We constructed a simple stochastic microscopic model of the growth of rigid multi- 
filament biopolymers. It was argued that association/dissociation rates of individual 
monomers depend strongly on local environment, which leads to the conclusion that there is 
a finite number of polymer configurations that specify the dynamics of growing biopolymers. 
We suggested that most relevant configurations have protofilaments at distances less than 
a monomer length from each other, i.e. "one-layer" configurations. As a result, the mean 
growth velocity and dispersion, or the effective diffusion constant, was calculated exactly for 
any number of protofilaments and for any shifts between them in terms of rate constants for 
attachments and detachments. It should be noted that simultaneous knowledge of the veloc- 
ity and dispersion provides a better description of fluctuations and variability in biopolymer 
dynamics. Precise experimental measurements of these properties would provide a valuable 
information on the growth mechanisms. 

Our approximate theoretical approach easily takes into account the effect of external 
forces on growing biopolymers. It could be done by modifying rate constants using chemical- 
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kinetic arguments. It allowed us then to construct a force-velocity relation and to estimate 
the stalling force, i.e., the force when the growth velocity becomes zero. It was suggested 
that the comparison with experimental force- velocity curves would provide a testing ground 
for this theoretical method. 

Explicit expressions for the mean growth velocity and for dispersion obtained in the "one- 
layer" model allowed us to investigate the dependence of growth processes on monomers 
concentration. It is found that at large concentrations the mean growth velocity grows 
linearly, in agreement with phenomenological descriptions. However, at low concentrations 
the significant deviations from linearity may be found for some sets of parameters. This 
observation contradicts the main result of phenomenological models which assume that 
linear dependence of the growth velocity is valid at all concentrations. These non-linear 
deviations are probably due to the fact that the "one-layer" model provides a more realistic 
microscopic description of structural and geometrical properties of growing biopolymers, 
which apparently is more important at low concentrations of monomers. This question 
requires careful experimental and theoretical tests. 

The vahdity and applicabihty of the "one-layer" model was discussed for the simple case 
of the growth of biopolymers with N = 2 protofilaments. For this case an exact solution 
for full dynamics, that accounts for all possible polymer configurations, was derived. It 
is found that the predictions of the approximate theory for the mean growth velocity are 
approaching the exact values for large, but realistic values of lateral interactions between the 
protofilaments. This indicates that the "one-layer" model probably captures main physical 
and chemical properties of complex growth processes, and it can be used to describe real 
systems, such as microtubules and actin filaments. 

The advantages of using the "one-layer" model to describe the growth of rigid biopoly- 
mers is not only its simplicity and the ability to obtain explicit expressions for dynamic 
parameters, but also its very fiexible and the fact that it can be extended and modified in 
several directions in order to describe these complex processes more realistically. Firstly, 
more polymer configurations may be included by considering a "two- layers" condition, i.e., 
that the protofilaments in relevant configurations are at distances less than two monomer 
subunit lengths. More layers can be added in a similar fashion if required. Thus the results 
can be improved iterativcly, i.e., this extension is analogous to a scries expansion approach 
to the full dynamics description. Secondly, the apparently weaker interactions at the poly- 
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mer lattice seam can be also incorporated. It will be interesting to know how this addition 
will effect the growth dynamics. Another possibihty is to include tubuhn-GTP hydrolysis in 
growth dynamics. This may give a new route to investigate the dynamic instability phenom- 
ena, which is still the most outstanding problem in the microtubule dynamics. Also, it will 
be interesting to compare the "one-layer" model with computer simulations of microtubule 
dynamics, which we plan to do in a future work. 
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Figure Captions: 

Fig. 1 A typical configuration of growing rigid biopolymer consisting of protofilaments. 
Tlie seam in the polymer lattice is between the protofilaments and 1. The tip of leading 
protofilament is at the origin. The free monomer can bind to any protofilament. 

Fig. 2 a) A biopolymer configuration with all protofilaments at distances less than one 
subunit length d from the leading protofilament, i.e., a "one-layer" configuration, b) The 
biopolymer configuration which is not a "one-layer" configuration. The protofilament 2 is 
at the distance 02 > d. 

Fig. 3 Four possible situations in the growth of a biopolymer with N = 2 protofilaments. 

Fig. 4 A kinetic scheme for full dynamics description of the growth of a biopolymer with 
N = 2 protofilaments. 

Fig. 5 Ratio of mean growth velocities as a function of the parameter 7 for different protofil- 
ament shifts. 

Fig. 6 The force-velocity curve from the experimental data of Dogterom and Yurke.— The 
solid line represents the fit obtained using the "one-layer" model with load-dependent rate 
constants. 
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